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ABSTRACT 



We report on near-infrared (IR) interferometric observations of the double- 
lined pre-main sequence (PMS) binary system DQ Tau. We model these data 
wi th a visual orb i t for DQ Tau supported by the spectroscopic orbit & analysis 
of iMathieu et all (119971 ). Further, DQ Tau exhibits significant near-IR excess; 
modeling our data requires inclusion of near-IR light from an 'excess' source. 
Remarkably the excess source is resolved in our data, similar in scale to the binary 
itself (~ 0.2 AU at apastron), rather than the larger cir c umbin ary disk (~ .4 



AU ra dius). Our observations support the IMathieu et all (119971 ) and ICarr et al. 



(120011 ) inference of significant warm material near the DQ Tau binary. 



Subject headings: binaries: spectroscopic — stars: circumstellar matter — stars: 
pre-main sequence — stars: individual (DQ Tau) 



1. Introduction 



DQ Tau (HBC 72) is among a small number of known "classical" T-Tauri (CTTS) 
spectroscopi c binaries. DQ Tau's strong Ha emission and de facto T Tauri status was 
reported by IJoyl (119491 ) . The system has a spectral ene rgy distribution ( S EP) typical of 
CTTS, with strong near and mid-infrared (IR) excess (jStrom et all Il989l : ISkrutskie et al 
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19901 ; iMathieu et al 19971. §12]). Ha emission and continuum veiling indicate s ig nificant accre- 



tion o nto the stars (IValenti et allll993l ; lHartigan et allll995l ; iBasri et al.lll997l ). IMathieu et al 
(119971 . herein M1997) established DQ Tau as an eccentric, short-period (15.8 d), double-lined 
spectroscopic binary composed of similar late K-stars; pres umably the mid-I R excess is due 
to emission from a circumbinary disk. Further, M1997 and lBasri et al.l (119971 . herein B1997) 
identified photometric and spectroscopic variability at the orbit period and phased near 
periastron, interpreting these variations as enhanced accretion as the stell ar components en- 
counter streams of infalling material from the circumbinary disk. Recently ISalter et all (120081 ) 
reported a mm flare in DQ Tau, but interpret the flare as interacting stellar magnetospheres 
similar to that seen in V773 Tau A (IMassi et all 120081 ). 



Binary systems are expected to clear inner gaps in their circumbinary disks (out to 
several times the binary semi-majo r axis) as material is d ynamically ejected by the compo- 
nents (lArtymowicz fc Lubowlll994l ; iPichardo et all 120051 ) . IJensen fc Mathieul (119971 ) report 
SED evidence for such inner-disk clearings in a set of T Tauri binaries. The SEDs in these 
systems show a deficit of near-IR emission relative to mid-IR flux - interpreted as a lack of 
the warmer inner-disk material that has been dispersed by t he stars' orbital motion. Con- 
versely, th e DQ Tau SEP (a nd a few others, e.g. AK Sco (I Jensen fc Mathieul 1 199 71 ). and 
UZ Tau E IJensen et all ( 120071 )) exhibits no such near-IR deficit, leading M1997 to conclude 
"there is clearly warm material within the binary orbit." M1997 modeled the DQ Tau 
SED with a modest (5x l0~ 10 Mm) amoun t of warm (1000 K), optically-thin material in the 
DQ Tau binary region. ICarr et al.l (120011 . herein C2001) supported this conclusion with IR 
spectroscopic detection of warm (1200 K) CO gas emission from DQ Tau. 

The presence of material inside the expected DQ Tau dynamical gap is unsurprising: the 
system shows strong accretion diagnostics. But the amount and morphology of this 'inner' 
material bears on how such PMS binary systems interact with and accrete material from their 
circumbinary reservoir, and motivates study with the highest-resolution techni ques available. 
Here we report on observations of DQ Tau with the Keck Interferometer (KI, IColavita et al 
20031 ). These observations partially resolve the DQ Tau system, and allow us to model the 
system visual orbit based on the spectroscopic orbit by M1997. Further, we find that we 
must account for a static 'excess' flux which is compact, but partially resolved in these data. 
We interpret this excess source with possible morphological models for the warm material 
postulated by M1997 and C2001, and discuss implications on inner material and accretion 
in the DQ Tau system. 
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2. DQ Tau Spectral Energy Distribution 



The DQ Tau S EP has been studied by IStrom et all (jl989|), ISkrutskie et all (Il990( l. 
Hartigan et all (119951 ). and M1997; as it bears on our KI data analysis we summarize the 
SED here. Figure [1] shows a visible and near-IR SEP mode l for the DQ T au stellar pho- 
tosph eres derived using photometry from IStrom et all (119891 ) and 2MASS (ISkrutskie et al 
2006 ). Adopting published extinction estimates (Ay = 2.13; Strom et al 19891 ) and photo- 
spheric parameters (T e g = 4000 K, log g = 4.0, solar abundance; M1997), we can model 
photometry between 0.45 and 1 /iin (excluding blue accretion flux and IR excess) with a 
single photosphere (the stellar components of DQ Tau are nearly identical; M1997). We find 
the net photosphere model shown in Figure [TJ corresponding to a total st ellar luminosity of 



0.88 L for D = 140 pc, in good agreement with previous estimates (e.g. IStrom et allll989 
M1997). However, these estimates ignore the significant (0.35 ± 0.15 at 0.6 /im) veiling 
reported by B1997. An alternative model assuming uniform veiling between 0.45 and 1 /im 
is also given in Figure [TJ corresponding to total stellar luminosity of 0.65 L . The total 
luminosity of the DQ Tau photospheres is likely between these two extremes. 

DQ Tau's SED exhibits strong IR excess beyond 1 /im. M1997 modeled this excess 
with a power-law, and particularly argued for the presence of a cooler circumbinary disk and 
additional warm material within the binary orbit based on substantial excess over the stellar 
emission between 1-5 /im. Most significant for our purposes is the excess in i^-band (2.2 
/xm) - highlighted in Fig [TJ Based on the photospheric models from above we estimate the 
fractional excess r at K over the stellar contribution to be in the range of 0.5 -1.1 (modulo 
intrinsic /^-variability of the system; we return to this que stio n in § HI) . The bottom of this 



range agrees with similar estimates from lStrom et all (119891 ) and ISkrutskie et all (119901 ) . while 
higher values result from encorporating the veiling estimate of B1997. This SED model makes 
it clear our KI observations will contain flux from both the DQ Tau binary components and 
additional material, but the exact value of the excess .fT-emission over stellar is uncertain to 
a factor of two. 



3. Observations and Orbital Modeling 

KI Observations The KI observable used for these measurements is the fringe contrast 
or visibility (specifically, power-normalized visibility modulus squared, V 2 ) of an observed 
brightness distribution on the sky. KI observed DQ Tau in K-band on five nights between 
2005 Oct 25 and 2007 Oct 28, a data set spanning roughly two years and 49 orbital peri- 
ods. DQ Tau and calibration objects were typically observed multiple times during each of 
these nights, and each observation (scan) was approximately 130 sec long. As in previous 
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Fig. 1.— DQ Tau SED Model. We model the DQ Tau photospheres with a 4000 K template 
(M1997) and significant extinction (Ay = 2.13; iStrom et all Il989l ). The red curve shows 
an alternative photospheric model assuming a veiling fraction of 0.35 ± 0.15 (B1997). The 
fractional TT-excess over the photosphere (r) ranges between 0.5 and 1.1 depending on which 
photosphere model is adopted. 
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public ations, our KI V 2 calibration follows standard procedures described in IColavita et al 
( 120031 ). For this analysis we use HD 27282 (G8 V) as our calibration object, resulting in 26 
calibrated visibility scans on DQ Tau over five epochs. The V 2 observations are depicted in 
Figure [2], along with our best-fit orbit model (discuss e d bel ow). Orbital analysis methods 
for such V 2 observations are discussed in lBoden et al.l (120001 ) and not repeated here. 



Figure [2] shows the DQ Tau calibrated visibilities are significantly less than one in all 
five epochs - indicating resolved structure in our observations. Further, the five epochs show 
variability with time and hour angle - indicating the source appearance changes with time 
(e.g. with phase of the DQ Tau binary orbit), and is non-axisymmetric (both as expected 
for a resolved binary system). 



Boden et al. 


2000, 


2005, 


2007) 



Orbit Model As in previous analyses (e.g. 
KI visibilities with a binary source. The DQ Tau physical parameters from M1997 and 
putative 140 pc distance imply an apparent semi-major axis on the order of 1 mas. This 
apparent separation is marginally resolved in our data (projected fringe spacing of 5.1 mas), 
and does not allow an independent solution for the binary visual orbit. Therefore we have 
co nstrained our orbita l modeling with parameters from M1997 (with period slightly revised 



by lHuerta et al.l 120051 ) . Effectively we solve only for Q and the sense of rotation on the sky 
(i.e. whether i is greater or less than 90°). Further, from § [2] it is necessary to account for 
the visibility contribution of DQ Tau's K excess (even if that flux is incoherent on angular 
scales measured in these data). We find it sufficient to model the visibility due to the excess 
flux with a single, static visibility offset parameter V£ Sset (discussed in §1]). 

Figure [3] depicts our relative visual orbit model, with the primary rendered at the origin, 
and the secondary rendered at the five orbit phases of our KI data. Figure [2] shows the 
comparison of the KI data and predictions; the model clearly matches the data well within 
the estimated errors. Table [T] summarizes our DQ Tau orbit model as adopted from M1997 
and derived here. Our data show that the orbit motion is clockwise on the sky (retrograde); 
we adopt an inclination value (157°) constrained by the M1997 sini estimate and reflecting 
retrograde orbit motion. 



4. Interpreting the Visibility Offset 

Offset accounts for contributions from DQ Tau's K excess in our modeling. We have 
taken this visibility offset as static, that is invariant in time and with baseline projection 
angle; this construct warrants some consideration. First we note that M1997 and B1997 re- 
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Fig. 2. — KI Data/Model Comparisons for DQ Tau V 2 Observations. Here we give com- 
parisons for our five epochs of KI V 2 observations and our best-fit orbital model with and 
without the estimated visibility offset. In each frame both KI data and model are shown. 
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Table 1. Orbital Parameters for DQ Tau. 



Orbital 
Parameter 


M1997 


This Work 


Period (d) 


15.8016 




T (MJD) 


49582.04 




e 


0.556 




K Aa (km s _1 ) 


21.6 




K Ab (km s _1 ) 


22.4 




7 (km s _1 ) 


22.4 




(deg) 


228 




Q (deg) 




179 ± 10 


z (deg) 


23 


157 


a (mas) 




0.96 


A K (mag) 







V 2 

v offset 




0.15 ± 0.03 



Note. Our orbit model is con- 

strained to parameters estimated 
by M1997 e xcept for the period 

V 2 

v offset ' 



flHuerta et al.l 120051 ) . fi and 



Note that the inclination i is constrained 
to the M1997 value, except that our data 
& modeling indicate the motion of the 
binary is clockwise (retrograde) on the 
plane of the sky. 
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Fig. 3. — Relative Visual Orbit of DQ Tau. The relative visual orbit of the DQ Tau system 
is given, with the primary rendered at the origin, and the secondary position shown at the 
five epochs/orbit phases of our observations (two of our epochs are at periastron). 



ported photometric variability in DQ Tau, but this was strongest in blue (U & B) colors and 
attributed to variability in accretion luminosity; as we will argue below spectrophotometric 
monitoring suggests changes in the system K flux are modest. Second, the assumption of an 
offset that is independent of baseline projection de facto assumes that the K excess centroid 
is centered (co-axial) on the binary center of mass/ligh t (for an equal- mass binary), and is 
itself axially symmetric. lArtymowicz fc Lubowl (119961 ) have modeled accretion in systems 
such as DQ Tau, and their modeling clearly shows time- variable and non-axisymmetric struc- 
tures. These time-variable features call into question our static visibility offset construct, 
but the degree of visibility variability will depend on the relative intensity of symmetric and 
non-symmetric components of the K flux in the DQ Tau excess, and the spatial frequency 
coverage in our interferometric data. In our last epoch (2007 Oct 28; Figure [2] bottom) we 
specifically made the observation at periastron, and pushed the KI instrument to its maxi- 
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mum range in hour angle (roughly four hours) to test this static offset / axisymmetic emission 
assumption to the greatest extent possible. Data from this last epoch show no significant 
signs of variation with hour angle, supporting the axisymmetric emission construct. Given 
these considerations we settle on a static visibility offset to account for the K excess because 
nothing more complicated is justified by our data: the binary plus static offset construct 
adequately model our KI data over their range of spatial frequency. 

The (monochromatic) V 2 of a two-component composite scene is given by: 

2 _ VI + r 2 K, 2 + 2rW 2 cos ^ {V x + rV 2 f 

composite ~~ M _|_ r \2 ~^ M _|_ r \2 

with V\ and V2 the visibilities of two components, r the flux ratio (2 to 1), and the phase 
difference between the two fringes (e.g. = 27r/AB-s for a typical binary source); the 
second form assumes this phase offset is zero or the two sources are co-axial. Identifying V\ 
and V2 with the DQ Tau binary and the i^-excess respectively, and evaluating the expression 
at a convenient binary phase (i.e. periastron, where Hinary ~ 1) leaves: 



t/2 / \ (1 "I" '"^excess) * t r2 

Composite (Periastron) « — — « 1 - V oSset 



This allows us to estimate the net visibility of the DQ Tau ^-excess in terms of V 2 ffset 
(Table OQ) as: 

^exeess = ^((l-V; 2 frset )( 1 / 2 )(l+r)-l) (1) 

Evaluating Eq. [1] with r in the range of 0.5 - 1.1 (§ [2]) and V^ Sset = 0.15 (Table [Q yields 
Kxeess = 0.77 - 0.85. 

M1997 documents significant optical variability in DQ Tau, and attributes the vari- 
ability to enhanced accretion near periastron. Pertinent to this discussion is the degree of 
K variability in th e system. DQ Tau has been monitored with the CorMASS instrument 



(Wi lson et alll200ll ) as part of a T Tauri accretion variability program discussed in lBary et al. 
(120081 ). Analysis of these data suggests that K emission during enhanced accretion is ap- 
proximately 20% brighter than during a quiescent phase. This variation is small compared 
to the factor of two uncertainty in r from the DQ Tau photospheric uncertainty. 



Remarkably this range for Kxcess indicates the KI data resolve, but do not over- resolve, 
the K excess, and this result is robust even incorporating the large uncertainty in r. For 
instance, if the dominant ii'-excess came from a warm inner-edge of the circumbinary disk 



(nom inally at a barycentric distance of ~ 0.4 AU: lArtymowicz &: Lubowlll994! : lPichardo et al 



20051 ). the apparent inner-edge diameter would be ~ 5.7 mas. The V of such a ring mor- 



phology would be ~ 0.1 - 0.2. Apparently the K excess is dominated by emission from a 
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region significantly more compact than the dynamically-allowed inner edge of the circumbi- 
nary disk. Our visibility measurement supports the M1997 and C2001 conclusion that there 
is warm material in the vicinity of the DQ Tau binary, and emission from that material 
dominates the system's K excess. 

If we neglect .fT-emission from the circumbinary disk entirely (the equilibrium temper- 
ature of material at 0.4 AU is on the order of 475 K), we can estimate the angular scale of 
emission from the warm material. Table [2] summarizes a small set of simple assumed emis- 
sion morphologies for the .fT-excess: Gaussian profile, thin ring, and uniform disk. Under 
these assumptions the characteristic apparent size required to match the Kxcess estimate is 
given. Of these three models we expect the Gaussian profile most closely matches the warm 
material emission profile; we note that this assumption results in a size scale (1.1 - 1.4 mas; 
0.15 - 0.2 AU at 140 pc) that closely matched the apparent binary separation at apastron 
^1.5 mas). Inc oherent scattered-light from the circumbinary disk (proposed for CTTS by 



Pinte et all 120081 ) would make this warm material characteristic size estimate smaller. 



5. Discussion 

We have modeled our DQ Tau observations based on orbital parameters from M1997 
and system SED. Unsurprisingly, the significant .ff -excess must be included in the V 2 mod- 
eling. Remarkably these data indicate the excess must come from a region smaller than the 
circumbinary disk. Further, our data suggest that this excess is distributed on the physical 
scale of the binary orbit (~ 0.1 - 0.2 AU) rather than being either much smaller (e.g. cir- 
cumstellar disks) or much larger (i.e. the circumbinary disk) than the stellar separation. 

Our data and modeling support the M1997 and C2001 inference that DQ Tau has sig- 
nificant warm material in the inner orbit region (in addition to the substantial circumbinary 
disk). Binary dynamics, accretion, and wind/outflow processes work to dissipate this inner 
material over a few binary orbital periods. The static visibility offset that adequately models 
our KI data over many DQ Tau orbital periods suggests the system is in quasi-equilibrium 
with material inflow replenishing dissipated material in the binary region. 

It is important to note that our present data on and modeling of this remarkable system 
is relatively crude. In particular the ~ 5 mas KI fringe spacing only partially resolves the 
~ 1 mas DQ Tau binary orbit. Clearly limited spatial information leads to limitations in 
our DQ Tau modeling, and some care in interpreting our conclusions is warranted. Our 
construct that the i^-excess morphology is axisymmetric see ms the most suspect; it run s 



counter to existing modeling of accreting binary systems (e.g. lArtymowicz fc Lubowlll996l ). 
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and photometric/accretion brightening near periastron (M1997). However nothing more 
sophisticated is justified by our data - in particular the KI data from our last epoch (2007 
Oct 28) test the axisymmetric modeling assumption to the practical limits of KI capabilities 
(in broad-band data). But there should be an asymmetric component to the near-IR flux at 
some contrast level and spatial scale in the DQ Tau system. Going forward it is important 
to probe the extent of any axisymmetry in the near-IR excess at a greater diversity of spatial 
scales and wavelengths to understand the distribution and flow of material in DQ Tau and 
similar systems. 

Data presented herein were obtained at the W.M. Keck Observatory. We gratefully 
acknowledge the support of personnel at the Jet Propulsion Laboratory, W.M. Keck Obser- 
vatory, and the NASA Exoplanet Science Center in KI observations of DQ Tau. We thank 
the anonymous referee for his comments on including the veiling luminoisty in the SED 
modeling for DQ Tau. 
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Table 2. Characteristic Angular Scales for Simple Morphological Models of DQ Tau 

if-band Excess. 



Model 


Characteristic 


Note 




Scale (mas) 




Gaussian 


1.1 - 1.4 


Gaussian FWHM 


Ring 


1.3 - 1.7 


Ring Diameter 


Uniform Disk 


1.9 - 2.4 


UD Diameter 



Note. — Assuming V; xcess = 0.85-0.77 
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